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Abstract 

describe  results  of  transient  two-dimension¬ 
al  finite  element  computations  used  to  simulate 
elastic  and  visco-elastic  surface  response  to  the 
stress  generated  by  the  surface  turbulent  boundary 
layer  in  water.  He  are  characterizing  the  tempor¬ 
al  and  spatial  response  of  soft,  pliable  surface 
coating  materials  such  as  natural  rubbers,  artifi¬ 
cially  treated  rubber-like  compounds,  and  cured 
polymeric  gels.  He  examine  a  variety  of  coatings 
from  stacked  homogeneous  layers  to  internally 
structured  layers  with  open  and  fluid  filled  cells 
or  with  directional  ribs  and  channels.  In  addi¬ 
tion,  we  examine  non-linear  visco-elastic  loss 
characteristics  and  time-dependent  influences  on 
material  response.  Our  goal  is  identifying  drag 
and  characterizing  the  most  important  dynamic 
material  properties  to  achieve  this  end.  Two 
dimensional  transient  numerical  computations  are 
applied  to  simulate  the  space  and  time  random 
scalar  pressures  which  excite  the  material 
response. 

IV 

'  I.  Introduction 

Our  investigation  focuses  on  compliant  coat¬ 
ings  which  exhibit  favorable  dynamic  response  rel¬ 
ative  to  critical  turbulent  boundary  layer  motions. 
Favorable  response  refers  to  surface  motions  which 
may  reduce  or  cancel  the  boundary  layer  motions 
and  consequently  reduce  skin  friction  drag.  The 
characterization  of  fluid-solid  interface  dynamics 
with  relation  to  material  properties  and  identifi¬ 
cation  of  the  fundamental  fluid/solid  dynamic 
coupling  mechanisms  with  a  potential  for  reducing 
skin  friction  drag  are  the  basic  purposes  of  this 
research.  Emphasis  is  placed  on  materials  and 
coatings  whose  compliant  material  properties  and 
response  characteristics  can  be  altered  to  respond 
to  the  turbulence  by  modifying  internal  structure, 
or  composite  material  stiffness,  or  dynamically 
altering  viscous  filler  properties  or  by  internal 
heating.  Such  coatings  may  be  effective  over  a 
broad  range  of  flow  conditions.  Our  analysis  is 
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specifically  directed  to  drag  reduction  on  coated 
surfaces  moving  in  water. 

Post  World  Har  II  attention  on  compliant  coat¬ 
ed  drag  reduction  possibilities  was  stimulated  by 
study  of  the  swimming  efficiency  of  mammals,  par¬ 
ticularly  porpoises.  One  of  the  early  post  World 
Har  II  investigations,  in  which  trial  coatings 
were  devised,  was  apparently  based  on  study  of  the 
porpoise  skin  surface  and  sub-surface  bone,  carti¬ 
lage,  fat,  and  vascular  composition  and  structure. 
These  coatings  were  developed  and  tested  in  the 
United  States  by  M.  O.  Kramer  in  the  1950's  and 
1960*3.1*2  Unfortunately,  the  experimental 
results  which  indicated  that  significant  drag 
reduction  could  be  achieved  using  a  combination  of 
viscoelastic  layered  materials  and  internal  struc¬ 
tures,  were  never  independently  and  reproducibly 
verified.  The  issue  remains  open. 

In  the  United  States,  subsequent  comprehensive 
studies  on  compliant  surface  drag  reduction  and 
materials  characterization  were  conducted  and 
reported  by  researchers  at  NASA  Langley  in  the 
1970*8.3***5  Suggestions  were  made  about  the 
mechanisms  by  which  transitional  boundary  layer 
stabilization  and  turbulent  boundary  layer  inten¬ 
sity  reduction  may  take  place.  The  notions  ad¬ 
vanced,  generally  focus  on  modulation  of  energetic, 
transient  boundary  layer  relatively  large  scale 
structural  features  which  have  been  observed  in 
comparatively  recent  turbulent  boundary  layer 
experiments.  Presumably  the  most  probable  drag 
reducing  dynamic  fluid  solid  coupling  mechanisms 
mav  be  suggested  by  these  experiments  and  further 
developed  by  theoretical  analysis  and  simulation 
prior  to  verification  by  carefully  defined  and 
conducted  experiments. 

Computer  simulations  of  the  complicated  bound¬ 
ary  layer  flow  and  instability  evolution  adjacent 
to  initially  prescribed  deformed 'surfaces  or  sur¬ 
face  motions  have  been  presented  by  Orszag.9*2®*l-l 
In  these  numerical  simulations  experimental  evi¬ 
dence  was  useful  in  suggesting  simplications  to 
the  computational  model.  By  modeling  certain, 
almost  periodic,  features  and  by  special  care  with 
respect  to  the  numerical  method  applied,  pseudo- 
spectral  computations  of  the  generalized  Navier- 
Stokes  model  solutions  are  obtained.  These  solu¬ 
tions  represent  the  most  sophisticated  turbulence 
simulations  adjacent  to  arbitrarily  deformed  sur¬ 
faces.  They  form  an  accurate  basis  for  analysis, 
provided  sufficient  attention  is  given  to  computa¬ 
tional  efficiency,  uniform  validity,  and  a  clearly 
defined  limit  to  application  range. 

He  make  frequent  use  of  the  results  and  inter¬ 
pretations  from  these  experimental  and  computa¬ 
tional  studies  in  our  computational  efforts  at 
identifying  potential  compliant  surface  materials 
and  drag  reducing  coatings.  The  most  promising 
are  projected  as  candidates  for  later  experimental 


study.  Some  of  our  preliminary  results ,  appearing 
here  for  background  information,  have  been  presen¬ 
ted  within  the  year.^ 

II.  Procedure 


bulent  boundary  layer  solutions  are  obtained  at 
the  start,  to  establish  the  necessary  estimated 
values  of  local  displacement  thickness,  dj' 
wall  shear,  £u,  root  mean  square  pressure 


To  develop  realistic,  experimentally  verified 
turbulent  excitation  pressure  space  and  time 
histories  we  make  use  of  one  of  our  co-author's 
(Ash)  two-dimensional  and  three-dimensional  Monte- 
Carlo  time  dependent  turbulent  pressure  simula¬ 
tions  to  generate  our  dynamic  loading  functions. 33 
Subsequently,  we  plan  to  explore  the  actual  coup¬ 
ling  mechanism(s)  using  the  general  numerical 
solutions  of  the  turbulent  Navier-Stokes  equations 
through  application  of  Orszag's  pseudo-spectral 
method. 9-11, 16  jn  thiS  paper  we  restrict  our 
discussion  to  an  intermediate  phase:  systematical 
materials  property  and  response  characterization. 
We  emphasize,  at  this  time,  thorough  evaluation  of 
compliant  materials,  special  internal  structures, 
and  passive  response  and  motion  modulation.  For 
this  we  apply  the  Ash  code  model  Monte  Carlo  pres¬ 
sure  simulations. £3 

We  consider  a  flat  plate,  with  zero  pressure 
gradient  immersed  immediately  below  the  surface  in 
sea  water,  we  examine  conditions  over  a  range  of 
specified  steady  state  free  stream  velocities. 

The  A3h  procedure!3  generates  a  modeled  scalar 
turbulent  boundary  randomly  fluctuating  surface 
pressure  using  a  quasi-Monte-Carlo  method  to 
develop  the  random  time  histories  at  preselected 
time  intervals  and  spatial  locations.  The  turbu¬ 
lent  pressure  development  is  considered  tripped  at 
a  selected  upstream  position  by  a  finite  (mm)  wire 
transition  trip.  A  hypothetical  development 
length,  50  wire  diameters  downstream  is  estimated 
for  initiation  of  the  fully  developed  turbulent 
boundary  layer.  This  defines  development  length 
and  point  of  origin  for  the  fully  developed  bound¬ 
ary  layer  fluctuating  pressure  simulations  genera¬ 
ted  by  the  Auh  code.  The  Ash  code  results  have 
been  shown  to  model  fluctuating  pressures  and 
scale  features  of  the  boundary  layer  in  close 
agreement  with  the  experimental  observations  of 
Bull14  and  of  Bull,  Thomas.35  Nevertheless, 
this  is  a  pressure  model  only.  Convection  and 
convection-dominated  bur3t-streak  and  upthrust 
vorticity  generation  features  are  separately 
imposed  using  observations  based  on  the  Offen- 
Kline  observations  associated  with  unsteady  char¬ 
acteristics  of  the  turbulent  boundary  layer.® 


<P'>  1  3-0*Cu  •  (3) 

And  time  scales  for  unsteady  burst  growth,  propo- 
gation  and  decay.  The  entire  burst  evolution  is 
assumed  convected  with  the  mean  flow  at  a  high  but 
constant  speed  proportional  .to  the  outer  edge 
velocity,  viz. 

uc  -  0.6  U.  (4) 

Experiments  conducted  at  the  NBS  and  Naval 
Research  Laboratories  on  promising  coatings  will 
verify  their  performance.  They  will  also  provide 
information  for  the  development  and  subsequent 
analysis  of  material  properties  and  the  crucial 
features  of  the  interactive  coupling  model. 

Next,  we  apply  a  discretized  procedure  to 
numerically  simulate  the  response  of  composite 
materials  and  structures.  To  accurately  follow 
high  frequency  motions,  in  a  practical  computation, 
both  implicit  and  explicit  two-dimensional  and 
three-dimensional  LLNL  finite  element  codes  of 
Hallquist  are  used. 17,18 

Among  the  principal  finite  element  code  fea¬ 
tures  consistently  U3ed  in  the  present  analysis, 
is  the  facility  for  treating  slide  lines  (or  in 
the  3®  version  sliding  surfaces)  between  unlike 
materials  or  at  geometric  discontinuities.  The 
interface  between  the  PVC  and  Neoprene  layers,  for 
example  is  treated  by  a  range  of  options  from 
friction  to  free  sliding  constraints  which  remove 
the  ambiguity  associated  with  discontinuous  mater¬ 
ial  properties  and  associated  erroneous  discontin¬ 
uities  in  the  continuum  description  of  the  dynamic 
stress  tensor  field. 

Generally,  the  unsteady  dynamical  problem  is 
formulated  in  a  lumped  mass  matrix  form  as  a 
perturbation  to  equilibrium  static  equations. 
Unsteady  dynamics  are  a  perturbation  to  the  equil¬ 
ibrium  system  which  is  relaxed  in  the  implicit 
NIKE  code  version  to  a  new  equilibrium  state  or 
time  marched  in  the  explicit  DYNA  code  version  to 
the  next  advanced  time  step.  The  lumped  mass 
matrix  form  is  written 


Specifically,  we  adopt  for  the  present  calcu¬ 
lation  the .Of fen-Kline  implications  of  the  depen¬ 
dence  of  the  spread  between  burst  events  as  pro¬ 
portional  to  the  inner  boundary  layer  (wall  region) 
variables. 

Ax  »  uQ  At  (1) 

wnere  Ax  i3  the  spatial  distance  between  burst 
events.  At  is  the  time  interval  and  u<j  is  the 
3hear  velocity  given  by  the  wall  shear,  £iu,  in 
accordance  with. 


Mdn+1  +  D(dn  +  aAd)  dn+1  +  K„(dn  +  oAd*)  Ad. 

-  P(dn)n+1  -  Fig")  (5) 

The  first  term  represents  the  nodel  mass  accel¬ 
eration  acting  on  a  consistent  lumped  mass  matrix, 
M.  Single  and  double  dots  over  the  deformation 
matrix  elements,  d,  refer  to  1st  and  2nd  deriva¬ 
tives  with  respect  to  time.  The  operator  matrix, 

£,  is  a  Rayleigh  damping  matrix. 


(2) 


D  =  a  ,M  +  Km 

as>  cJas  d  *T 


(6) 


where  p  is  the  fluid  density  and  y  is  the  normal 
coordinate  displacement  from  the  wall  surface 
through  the  boundary  layer.  Incompressible,  tur- 


where  £7,  is  the  material  property  dependent 
symmetric  positive  definite  stiffness  matrix,  J> 
represents  the  external  loads  vector  and  T  is  the 
stress  divergence  vector  centered  at  time  "n". 


/ 


The  constants  aj  and  8,j  are  usee  defined  at  initi¬ 
ation  of  a  specific  problem.  The  term  d*  repre¬ 
sents  the  estimated  displacement  vector  associated 
with  elemental  loads  and  acceleration  at  the  start 
of  an  iteration.  In  addition  to  the  slide  line 
and  slide  surface  features,  the  freedom  to  apply 
body  force  or  concentrated  nodal  loads,  pressure 
boundary  conditions,  shear  boundary  conditions, 
displacement  boundary  conditions  all  within  a  self 
consistent  Jaumann  stress  rate  formulation,  serve 
useful  purposes  in  our  applications.  In  our  simu¬ 
lations  we  are  specifically  concerned  with  compos¬ 
ite  structures  of  arbitrary  scale  size  in  dynamic 
simulations  with  problem  defined  time  constants. 

The  material  models  available  include  elastic, 
orthotropic  elastic,  elastic-plastic  and  linear 
visco-plastic.  The  first  and  last  two  are  of 
particular  interest  to  us  in  our  application. 

These  permit  modeling  of  variable  composite 
material  properties  and  constitutive  relations, 
and  modeling  of  the  internally  structured  compli¬ 
ant  surface  coatings,  such  as  those  with  air  gaps, 
water  passages,  and  rigid  reinforcing  stubs. 

These  interior  structures,  in  some  cases,  are 
found  to  effectively  shorten  the  characteristic 
response  wavelength.  Addition  of  heat  to  “tune* 
the  material  elastic  moduli,  and  specified  inter¬ 
nal  fluid  pressures  associated  with  internal  cell- 
to-cell  liquid  flow  (which  modifies  both  structur¬ 
al  stiffness  and  damping)  are  also  modeled,  quite 
generally.  The  present  finite  element  results  are 
for  pure  elastic  and  damped  viscoelastic  modeled 
composite  materials  in  transient,  plane  two-dimen¬ 
sional  strain.  The  three-dimensional  surface  is 
considered  composed  of  streamwise  infinitesimally 
thin  strips  excited  by  the  two  dimensional  trans¬ 
ient  random  pressure  oscillations.  Me  examine  the 
displacement  response  of  various  surface  materials, 
composites  and  structures.  Boundary  conditions 
include  both  clamped  and  roller  end  conditions  on 
the  strip  with  a  rigid,  non-yielding  lower  sur¬ 
face.  Full  three  dimensional  driven  surfaces  ace 
simulated  in  later  calculations,  together  with 
numerical  simulations  of  the  complete  Naviec- 
Stokes  equations  of  motion  following  the  pseudo- 
spectral  Orszag  procedure.  This  work  in  current 
progress  is  directed  to  definition  and  evaluation 
of  the  full  f luid-to-solid. to-fluid  coupling  and 
will  be  described  in  later  publications. 

III.  Results 

We  present  preliminary  results  of  our  previous 
computations  and  analysis'-^  to  assist  descrip¬ 
tion  of  our  numerical  procedures  and  the  physical 
interpretation  of  our  currently  available  results 
with  respect  to  potential  drag  reduction. 

Figures  1  and  2  show  elastic  response  of  indi¬ 
cated  materials  as  characterized  by  response  amp¬ 
litude  and  circular  frequency,  respectively.  Both 
sets  of  calculations  were  made  for  zero  immersion 
depth  ocean  surface  (Z  »  0)  motion.  The  solid 
curves  are  the  response  envelopes  for  the  indica¬ 
ted  speeds:  10,  20,  and  10  m/s.  The  dashed 
curves  intersecting  ace  the  response  intercepts 
for  the  specified  materials:  natural  rubber,  neo¬ 
prene  and  plastisol  gel  (PVC) .  The  amplitude  res¬ 
ponse,  d.  Fig.  1  is  ratioed  to  a  critical  boundary 
layer  dimension,  Sj'  (displacement  thickness)  rel¬ 
ative  to  dimensions  of  potentially  unstable  motions 
of  the  wall.  The  rms  turbulent  presture  fluctua¬ 


tion  <p‘>  is  ratioed  to  the  Young's  modulus  (E)  of 
the  materials  in  these  illustrations. 

Figure  3  illustrates  configuration,  computa¬ 
tional  geometry,  and  displacement  response  results 
for  elastic  analysis  of  homogeneous  layers  of  PVC 
on  neoprene  and  the  reverse  ordering  of  the  layers, 
neoprene  on  PVC.  Pure  sinusoidal  excitation  pres¬ 
sures  were  used  for  these  preliminary  studies.  We 
found  that  in  this  simple  elastic  limit  the  dis¬ 
placements,  d,  ratioed  to  the  depth  of  material 
layer,  h,  is  functionally  characterized  by  a 
dimensionless  group  relationship, 

£  -  constant  *  (l-v2)  (^)  f(-— )  .  (7) 

Here  v  is  the  composite,  effective  Poisson's 
ratio,  b  is  the  effective  pressure  patch  breadth, 

<i>0  is  the  characteristic  natural  frequency  of  the 
material  layer.  The  layers  show  the  normal  tend¬ 
ency  to  amplitude  divergence  indicated  near  the 
Characteristic  frequency,  uig.  The  softer  PVC  on 
the  upper  surface  is  clearly  superior  to  the  oppo¬ 
site  arrangement  for  obtaining  the  largest  ampli¬ 
tude  response  in  material  layers.  Similar  dimension¬ 
al  analysis  collapse  of  the  response  histories  for 
surface  coatings  consisting  of  visco-elastic 
elements  has,  so  far,  been  unsuccessful.  However, 
such  a  characterization  is  sought  from  our  analysis 
in  order  to  provide  a  basis  for  material  and 
structural  characterization  when  subjected  to  ran¬ 
dom  turbulent  loads. 

The  non-linear  frequency  dependence  of  the 
shear  modulus  for  one  of  the  materials,  "plastisol" 
(PVC)  is  shown  in  Fig.  4.  This  experimental 
(dashed  curve)  material  frequency  dependence  is 
modeled  satisfactorily  (solid  curve)  in  the  finite 
element  plane  strain  calculations.  The  experimen¬ 
tal  data  was  provided  to  us  by  researchers  at  the 
Naval  Research  Laboratory.  Figure  5  shows  the 
relative  variation  in  amplitude  response  for  the 
viscoelastic,  Internally  structured  material  coat¬ 
ings  possessing  internal  cells  of  about  1/3  total 
coating  depth.  They  are  numerically  tested  with 
various  "fillers"  in  the  open  channels.  Fillers 
included  tests  with  void  (actually  air  at  1  atm.) , 
plastisol  filler-  water,  and  neoprene  filler  with 
the  order  of  response  magnitude  following  the  same 
order  as  the  listed  fillers.  Voids  give  the 
greatest  displacement  and  neoprene  filler  the 
least.  However,  considering  both  response 
frequency  and  amplitude,  the  visco  elastic  incom¬ 
pressible  fillers  are  the  most  promising  of  the 
materials  simulated  in  these  numerical  tests. 

Three  classes  of  materials,  geometries  and 
structures  for  compliant  coatings  are  studied  here. 
We  examine,  in  turn,  homogeneous  multiple  material 
layered  coatings,  internally  structured  coatings 
with  imbedded  cells  or  channels  filled  with  visco¬ 
elastic  incompressible  fluid,  and  a  small  dimen¬ 
sion  (2  mm  thick)  Kramer  coating  suggested  by  his 
later  studies. ^ 

Figures  6  through  12  sunnarize  our  findings  on 
uniformly  homogeneous  layered  structures  with  10% 
PVC  plastisol  overlaying  a  neoprene  sub-surface 
layer  in  a  compliant  coating  over  a  relatively 
rigid  (aluminum)  surface.  The  geometry  and  sur¬ 
face  deformation  pattern  (enlarged  by  the  explicit 
scale  factor)  is  shown  at  time  zero  (static)  and 


z 


thretf  subsequent  tines,  where  times  are  given  in 
sec. 


The  driving  pressure  <Pa)  tine  histories  at 
25%,  50%  and  75%  stations  from  the  leading  edge  of 
the  compliant  surface  are  shown  in  Fig.  7.  The 
corresponding  Fourier  transformed  pressure  spectra 
are  shown  in  Fig.  S.  Peak  frequency  shifts  to 
lower  values  as  the  wave  train  energy  is  distribu¬ 
ted  over  a  larger  streaawise  dimension  and  decays 
with  increasing  involvement  in  the  modeled  bound¬ 
ary  layer  fluid. 

The  resulting  surface  displacement  amplitude 
time  histories  are  shown  in  Figs.  9  and  10  for 
free  stream  speeds  of  10  and  30  m/s,  respectively. 
The  amplitude  spectra  are  illustrated  in  Figs.  11 
and  12  for  the  10  and  10  m/s  speeds,  respectively. 
Shifts  of  the  peak  amplitude  spectra  to  lower 
frequency  than  that  for  the  driving  fluid  is  more 
pronounced  due  to  larger  solid  material  inertia, 
rigidity,  and  dispersive  characteristics.  From 
our  previous  elastic  and  viscoelastic  studies  the 
softest  (PVC)  materials  tested  is  usually  placed 
on  the  surface  for  developing  greater  response 
amplitudes. 

We  simulated  both  softer  and  harder  surface 
materials  over  structured  internal  layers  in  coat¬ 
ings  depicted  in  Figs.  13  through  22.  Figures  13 
and  14  show  geometry  of  the  structured  layers  and 
deformation  of  the  surface  at  designated  times 
(sec)  for  neoprene  and  PVC  surfaces  respectively. 
The  flow  speed  is  30  m/s.  Turbulent  pressure  time 
histories  and  pressure  spectra  are  shown  for  three 
designated  stations  in  percentage  of  length  from 
leading  edge  in  Figs.  15  and  16. 

The  PVC  on  the  upper  (fluid-solid  contact) 
surface  gives  the  largest  (most  promising)  dis¬ 
placement  amplitudes.  However,  the  frequency  of 
surface  motions  is  somewhat  lower  for  the  softer 
materials,  and  is  probably  too  low  for  useful  in¬ 
terference  with  the  oscillatory  boundary  layer 
structural  features.  The  displacement  amplitudes 
are  shown  for  10  m/s  and  30  m/s  in  Figs.  17  and 
18,  respectively  while  the  corresponding  amplitude 
spectra  are  shorn  in  Figs.  19  and  20.  The  shift 
of  the  peak  of  the  pressure  spectra  is  evidently  a 
function  of  increasing  length  of  coating  as  illus¬ 
trated  by  the  values  at  the  25%,  50%  and  75%  sta¬ 
tions  from  the  leading  edge  of  the  plate  model. 

Also  evident  is  the  lower  frequency  for  peak 
amplitude  of  response  spectra  in  comparison  to  the 
driving  spectra  shown  in  Fig.  16.  Here  again 
material  Inertia,  stiffness  and  dispersive  charac¬ 
teristics  are  dominating  influences  on  the  fluid 
to  solid  coupling  dynamics.  Significantly  more 
variation  in  response  will  be  possible  by  modula¬ 
tion  and  tuning  of  the  coating  structure  dimensions 
and  filler  characteristics.  These  variations  or 
permutations  are  more  conveniently  conducted  on 
the  computer  in  order  to  usefully  limit  and  spec¬ 
ialize  on  a  few  well-defined  candidates  for  sub¬ 
sequent  experiments. 

Figures  21  and  22  show  the  response  spectra  at 
10  and  30  m/s,  respectively  for  the  same  Internally 
structured  coatings  but  with  the  harder  neoprene 
replacing  the  softer  PVC  for  the  surface  layer 
with  the  PVC  moved  to  the  substrate  coating  layer 
position.  Amplitudes  are  significantly  reduced 
with  the  firmer  surface  materials.  However,  the 
slight  shift  to  higher  frequency  for  peak  response 


spectra  is  partially  modulated  by  the  characteris¬ 
tics  of  the  Internal  structure  and  sub-layer 
mater ial. 

Figures  23  through  26  show  the  latest  small 
scale  (higher  frequency  response)  Kramer  ribbed 
structures  suggested  by  his  later  work. 19  These 
simulations  were  made  at  a  lower  velocity  (1.5  m/s) 
suggested  by  his  experiments.  While  the  higher 
frequencies  are  excited  the  response  amplitudes 
(micron  or  less)  are  negligible  for  significant 
boundary  layer  interference  purposes.  Figure  23 
illustrates  the  Kramer  geometry  and  deformations 
at  selected  times.  Figure  24  shows  the  driving 
turbulent  pressure  spectra  at  1.5  m/s.  Figures  25 
and  26  show  this  Kramer  structure  response  ampli¬ 
tude  and  response  spectra,  respectively. 

IV.  Preliminary  Interpretations 

In  general,  the  response  of  all  mono  or  sequen¬ 
tially  layered  coatings  in  water  parallel  the 
suggestions  derived  from  the  gas  flow  predictions 
of  Bushnell,  Hefner  and  Ash. 5  Specifically 
unstructured  (non-cellular)  passive  (non-heated, 
non-driven  coatings)  unless  very  thick  (two  or 
three  times  substrate  thickness)  exhibited  little 
or  no  amplification  modulation  of  the  response 
from  the  driving  forcing  function  that  was  dis¬ 
tinct  from  that  of  the  rigid  metal  substrate. 

This  reponse  dominance  by  the  substrate  character¬ 
istics,  with  negligible  coating  modulation,  be¬ 
comes  more  pronounced  as  speed,  burst  frequency, 
and  simulated  depth  increased.  The  preliminary 
results  imply  that  for  phase  resonance  between 
forcing  oscillation  and  the  coating,  the  sub¬ 
membrane  cavity  would  have  to  be  filled  with  an 
incompressible  liquid  to  prevent  deformation 
extensive  enough  to  completely  overwhelm  other 
desirably  broad  band  structural  response  charac¬ 
teristics. 

Inertial  contrasts  between  fluid  and  solid 
coating  as  well  as  evident  solid  stiffness  signif¬ 
icantly  reduce  the  ratio  of  fluid  excitation  amp¬ 
litude  to  solid  response  amplitude  and  increase 
the  phase  lag  between  fluid  excitation  and  solid 
response.  These  factors,  in  conjunction  with 
natural  material  dispersion,  suggest  that  attain¬ 
ment  of  simple  fluid/solid  direct  interference 
coupling  is  improbable.  It  further  suggests  that 
for  potentially  effective  drag  reduction,  combina¬ 
tions  of  structure  and  materials  must  be  sought 
that  provide  multi-modal  response  to  pulsed  exci¬ 
tation  from  the  turbulent  boundary  layer.  One 
seeks,  in  this  way,  to  enlarge  the  probability  of 
generating  appropriately  phase  canceling  or  atten¬ 
uating  surface  motions. 
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Fig.  1  Elastic  theory  response  amplitudes  for 
selected  coating  materials  as  a  function  of  rms 
pressure. 


E 

Fig.  2  Elastic  theory  response  frequencies  for 
selected  coating  materials  as  a  function  of  rms 
pressure. 
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Fig.  8  Turbulent  pressure  spectra  at  10  m/s 
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Fig.  15  Turbulent  pre'ssure  at  30  m/s 


Fig.  13  Internally  structured  neoprene  surface 
coating  deformation  at  30  m/s. 
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Fig.  16  Turbulent  pressure  spectra  at  30  m/s 
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Fig.  14  Internally  structured  PVC  surface  coating 
deformations  at  30  m/s. 
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Fig.  19  Internally  structured  PVC  surface  coating 
Fig.  17  Internally  structured  PVC  surface  coating  response  spectra  at  10  m/s. 

displacement  response  at  10  m/s. 
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Fig.  18  Internally  structured  PVC  surface  coating 
displacement  response  at  30  m/s. 
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Fig.  20  Internally  structured  PVC  surface  coating 
response  spectra  at  30  m/s. 
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Fig.  21  Internally  structured  neoprene  surface 
coating  response  spectra  at  10  m/s. 
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Fig.  23  Internally  structured  Kramer  model  small 
scale  coating  and  modeled  surface  deformation  at 
1.5  m/s. 
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Fig.  24  Turbulent  pressure  spectra  at  1.5  m/s. 


Fig.  22  Internally  structured  neoprene  surface 
coating  response  spectra  at  30  m/s. 


Fig.  25  Internally  structured  Kramer  model  small 
scale  coating  displacement  response  at  1.5  ra/s. 


Fig.  26  Internally  structured  Kramer  model  small 
scale  coating  response  spectra  at  1.5  m/s. 
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